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Abstract 

We find that perturbation theory fails to describe the short range static potential 
as obtained from quenched lattice simulations, at least for source separations r > 
(6 GeV) -1 . The difference between the non-perturbatively determined potential and 
perturbation theory at short distance is well parameterised by a linear term with a 
slope of approximately 1.2 GeV 2 , that is significantly bigger than the string tension, 
a « 0.21 GeV 2 . 



1 Introduction 



The operator product expansion (OPE) constitutes the standard framework 
of including non-perturbative contributions to a QCD observable. Within the 
OPE, non-perturbative power corrections are included in the form of infrared 
sensitive operators of a given dimension that are accompanied by ultraviolet 
coefficient functions. These Wilson coefficients are in general themselves only 
fixed up to power corrections that can be reshuffled into higher dimensional 
operators without affecting physical results. This ambiguity is related to the 
presence of infrared renormalons [1] and also reflects the asymptotic character 
of the QCD perturbative series. 
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The powers of the correction terms within the OPE are completely determined 
by the symmetries of the operator that is expanded. In a standard infrared 
renormalon analysis the dimension of the lowest order renormalon contribution 
coincides with that of the lowest non-perturbative condensate that contributes. 
In the case of the running coupling this would for instance be the gluon con- 
densate [2], such that the leading order power correction to a s (q) is expected 
to be proportional to A 4 /g 4 . This picture has recently been challenged by the 
ultraviolet renormalon technique [3] and dispersion approaches [4-6]. It has 
been argued that in addition to the terms expected from the OPE, a further 
ultraviolet renormalon contribution, proportional to A 2 /q 2 might be present. 
Insofar, renormalons seem to allow for a more universal classification of power- 
like corrections than the OPE [7]. Power corrections to a s have been known 
to arise naturally in many physical schemes for quite a while [8]. 

Evidence of an unexpected A 2 /q 2 power correction to the gluon condensate has 
been obtained in the lattice study of Ref . [9] . Recently, the running coupling 
from the three-gluon vertex determined in Landau gauge on the lattice has 
been investigated and also here evidence for a contribution proportional to 
A 2 /q 2 has been reported [10]. 

In case of the static potential in position space, such a A 2 /q 2 contribution to 
the running coupling results in a term proportional to the quark separation, r, 
while the leading order correction from the standard renormalon analysis de- 
pends only quadratically on r [11]. Indeed, based on various reasonable model 
assumptions such linear short-range contributions have been obtained [11-13]. 

The present letter is motivated by these latest developments and the phe- 
nomenological relevance of the short-range static QCD potential in view of 
calculations of the it threshold cross-section. We summarise recent perturba- 
tive results, describe the method that we use to subtract the QCD perturbative 
series and, indeed, find the non-perturbative contribution to be dominated by 
a linear term. 



2 The potential in perturbation theory 



Recently, the static potential, 

V(g) = (1) 

has been calculated to two loops by Peter [14,15]. This calculation has been 
independently confirmed by Schroder [16] who corrected a normalisation error 
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of the original reference. The result reads, 

a v(q) = (*Ms(q) (l + a^Msiq) + a 2 a\^) 



(2) 



For the quenched approximation (ny = 0) they obtain, 



C(3) 1.202056903 denotes the Riemann (^-function, ai has already been 
calculated some time ago [17,18]. 

A Fourier transformation yields, 

V(r) = - Cf ^±, (4) 



for the potential in position space with [18,15,19], 

an(l/r) = a v (/j,) [1 + —^-ay) , ji = exp(-7 B )/r. (5) 



'-/e denotes the Euler constant. The coefficients of the QCD /^-function for 
iif = are given by [20], 

B- U 3 ~ 102 3 m - 2857 (6) 



By employing a recursive lattice finite size technique [21], the ALPHA collab- 
oration has recently obtained a value for the running coupling in quenched 
QCD [22], 

A^ = 0.602(48) t-0 1 . (7) 



r ~ 0.5 fm denotes the Sommer [23] scale of eq. (10) below. 

We choose to parametrise the solution of the renormalisation group equation 

by, 



I , , . 2 lnt d 2 dl\n 2 t A&lnt 



i + rfllni + < _ + __^_ + ft5 (2 



-i 



, (8) 
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with t = 21n(/x/Ajgg), d\ — |? and d 2 = dj — ||, which is consistent to 

order ccf. The integration constant corresponds to the standard normalisation 
convention of the QCD A-parameter. 



3 The lattice potential 



We have collected lattice data on the static potentials from Refs. [24-27]. All 
these results have been obtained in the quenched approximation with Wilson 
action at values of the bare lattice coupling 6.0 < (3 = Q/g^ < 6.8 that 
correspond to lattice spacings 0.19r > a > 0.069 r ~ 0.035 fm. The linear 
lattice extents are kept at a similar size in physical units and range from L = 16 
at P — 6.0 to L = 48 at (3 = 6.8. In all cases, we find La > 3r xs 1.5 fm, such 
that finite size effects are negligible [24,28]. 

All lattice data within this /5-range can be well fitted to the Cornell ansatz, 
aV{r)=aV {a)- -^-r + {aa 2 ){ra- 1 ) (9) 



for r > 0.4 r , where rotational invariance is found to be restored within 
statistical errors. The parameter V (a) that contains the self energy of the two 
static sources diverges like — l/[aln(Aa)] as the continuum limit is taken. The 
scale r is defined by [23], 



dV(r) 
dr 



1.65. 



(10) 



r=ro 



From quarkonium phenomenology one finds [23,26], r ~ 0.5 fm. We obtain 
this scale from our fits to eq. (9): r a _1 = \J (1.65 — e)/(cra 2 ). 

In Fig. 1, we plot the short range part of the lattice potentials obtained at 
different lattice spacings ap in units of r . The statistical errors are all smaller 
than the corresponding data symbols. The normalisation convention adopted 
to cancel the self-energy is V(r ) = 0. The first four data points measured at 
each /5-value correspond to r jap = 1, \/2, \/3, 2, respectively. The data for r > 
2a/3 all follow a universal continuous curve and rotational invariance is restored 
qualitatively. We take this as indication that for these distances, the continuum 
limit is effectively realised. For comparison, the result from the Cornell fit to 
r > 0.4 r , used to determine the scale r , is plotted (dashed curve). Towards 
small r the data points lie above this curve, indicating asymptotic freedom, 
i.e. a weakening of the effective Coulomb coupling with the distance. 
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Fig. 1. The quenched lattice potential versus perturbation theory in units of 
To ~ 0.5 fm. The only free parameter is the self-energy of the static sources. The 
lattice data has been normalised such that V(ro) = 0. The values 0.94 and 
0.77 Tq 1 have been subtraced from the perturbative one- and two-loop formulae, 
respectively, to allow for better comparison. 

In addition, we have included the one-loop (dashed dotted curve) and two-loop 
(solid curve with error band) continuum expectation of eqs. (2)-(8). Note that, 
up to the additive self-energy constant, no free parameter can be adjusted. 
For better comparison with the lattice data, normalised to V(r ) = 0, we 
have subtracted the values 0.94 Tq 1 and 0.77r _1 from the one- and two- loop 
expectations, respectively. The perturbative expression has a pole at A^- f» 
1.66 r ps (240 MeV) , such that we would not expect perturbation theory 
to reliably describe the lattice data for r > 0.5 r ps 0.25 fm or so. As can be 
seen from the plot, even at the smallest distance at our disposal, r ps 0.07 r ~ 
(24Ajg^-) _1 , the lattice data substantially deviates from perturbation theory; 
at any given distance the perturbative formula considerably under-estimates 
the slope dV/dr. 

In Ref. [29] it has been pointed out that due to the different form of higher 
order corrections, depending on whether one calculates the potential in po- 
sition or momentum space, numerically significant differences arise. We have 
convinced ourselves that in the quenched theory for 0.05 r < r < 0.4 r this 
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systematic uncertainty can indeed be numerically up to two times as big as the 
error due to the uncertainty in displayed in the plot. However, we found 
that the observed disagreement with perturbation theory cannot be repaired 
by fiddling around with the parametrisation. We remark that the difference 
between one- and two-loop perturbation theory is quite big as indicated by 
the value, $ = f3^ s + tuft + (a 2 - a?) A) ~ 3.19 > ^ s « 0.72. 

In position space, the correction is even more significant: starting from eqs. (2) 
and (5) and the definition, 

a R (l/r) = <*Ms(l/r) \l + da^l/r) + c 2 a\J , (11) 



we obtain, 

c 1 = a 1 + 2 lE (3 , c 2 = a 2 + 4 lE {3 oai + 2 lE (3 1 + ( 4 7 | + '— ) ffi. (12) 



This results in the numerical value, (3% = /i^ 5 + c 1 fi 1 + (c 2 — cj)/3 a* 6.70. 



4 Non-perturbative ultraviolet contributions 



Given the asymptotic character of QCD perturbative series, it is always am- 
biguous how to subtract perturbative contributions to a given physical process 
in order to isolate non-perturbative terms. Bearing this in mind, we employ 
the procedure of subtracting the one- and two-loop result, given by eqs. (2)- 
(8). Needless to say, that using a different scheme or cutting off the series at a 
different order in a s will in general yield somewhat different results. We would 
like to remark that beyond the two-loop level there appear to be difficulties 
in a consistent perturbative definition of the static potential [30,31]. 

In Figs. 2 and 3, we have subtracted the one- and two-loop continuum per- 
turbation theory result that correspond to the central value A-^ = 0.602 Tq 1 
from our lattice data. The resulting points are compared to linear curves (solid 
lines) with slopes 4.6 r^ 2 0.72 GeV 2 and 7.7 r^ 2 as 1.20 GeV 2 , respectively. 
We have added the error band due to the uncertainty in Ajjg to these curves. 
For r > 0.4 r and r > 0.3 r for the one- and two-loop results, respectively, 
deviations from the linear behaviour become visible while for r values that are 
small in units of the lattice spacings ap the points scatter around the inter- 
polating curve due to lattice artefacts. The data indicate the dominant short 
range correction to perturbation theory to be proportional to r. However, 
additional quadratic contributions might very well be present too. 
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Fig. 2. The quenched lattice potential with one-loop perturbation theory subtracted, 
in comparison with a linear curve with slope 4.6 r^ 2 . 



Given the visible lattice structure at small r, we do not attempt a fit to 
the data but we conservatively estimate systematic uncertainties of up to 
30 % on the observed slopes by allowing the interpolating curves to touch the 
most extreme combinations of data points that are possible. We checked that 
this error comfortably accomodates possible effects of additional quadratic 
contributions. The results are 3.4 ± 1.0 and 5.7 ± 1.7 times bigger than the 
string tension a m 1.35 Tq 2 for the one- and two-loop analysis, respectively. 
Neither a term of the size of the string tension as suggested in Ref. [32] nor a 
linear short-range term AV 2 (r) = ^a s ar predicted by Simonov [12] can alone 
account for this slope. Thus, additional non-perturbative linear contributions 
are likely to be present. 

We would like to emphasize that the independent determination of the QCD 
A-parameter by the ALPHA collaboration [22] was essential for our study: it 
is always possible to effectively compensate part of the slope by increasing the 
value of the A-parameter that enters the perturbative formula. 
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Fig. 3. The quenched lattice potential with two-loop perturbation theory subtracted, 
in comparison with a linear curve with slope 7.7 r^ 2 . 

We choose to parameterise the running coupling into, 

an(l/r) = a R ^-i 00 p(l/r) + c R r 2 (13) 



and find c R = —(0.90 ± 0.27) GeV 2 . A transformation into momentum space 
yields, ay(g) = av,2-ioo P (q) + cy/q 2 with cy = — cr/2. It is not a priori clear 
how to convert the result into different schemes and how to apply it to different 
physical observables. In fact, even the conversion to momentum space is not 
exact since the Fourier transformation affects higher order contributions and 
thus the difference between our result and perturbation theory. This can then 
result in a somewhat different slope of the linear term, due to the asymptotic 
character of the perturbative series. 
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5 Conclusion 



We find that perturbation theory fails to describe the short range static po- 
tential as obtained from quenched lattice simulations, at least for source sepa- 
rations bigger than (6 GeV) -1 . The difference between the non-perturbatively 
determined potential and perturbation theory at short distance is well parame- 
terised by a linear term. We quote the two-loop slope, (1.20±0.36) GeV 2 as our 
final result, which is significantly bigger than the string tension a 0.21 GeV 2 . 

We remark that the slope of the linear short distance contribution to the po- 
tential is, due to the large size of the two-loop coefficient, somewhat affected 
by the details of the resummation of the perturbative series. The precise value 
also depends on whether one allows for possible subleading non-perturbative 
contributions that are quadratic in r. This latter source of uncertainty is in- 
cluded into our error estimate. 

It is certainly worthwhile to repeat the present calculation using lattice actions 
different from the Wilson action to confirm universality of the result. A further 
significant reduction of the lattice spacing is at present ruled out by the size 
of available computers similar study with inclusion of sea quarks. In 

order to keep the argumentation transparent, we did not attempt to correct 
the potential for perturbatively known tree-level or one-loop lattice artefacts. 
This can in principle be done to reduce the amount of violations of rotational 
invariance in the short distance lattice potential. Unfortunately, no two-loop 
lattice perturbation theory result on the potential exists so far, such that a 
direct study of renormalon contributions on the lattice rather than in the 
continuum is not yet consistently possible to order a 3 . 

It should be interesting to investigate the same question in SU(N) gauge 
theory in three space-time dimensions (which is super-renormalisable) as well 
as to study the behaviour on small lattice volumes, where it is possible to 
further reduce the lattice spacing. As a consequence of an altered vacuum 
structure, confinement is lost on spatial volumes smaller than about 1 fm 3 . 
However, contributions that are related to ultraviolet physics should still be 
present and could easily be disentangled from terms that arise due to infrared 
properties of the theory. 
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